d-f Coulomb and quadrupole-strain interactions in DyB2C2 observed by resonant 

x-ray scattering 



T. Matsumura,^'0 D. Okuyama/ N. Oumi,^ K. Hirota,^ H. Nakao,^ Y. Murakami,^ and Y. Wakabayashi^ 

'Department of Physics, Graduate School of Science, Tohoku University, Sendat, 980-8578, Japan 
^Institute for Solid State Physics, The University of Tokyo, Kashiwanoha, Kashiwa, 277-8581, Japan 
^Institute of Materials Structure Science, High Energy Accelerator Research Organization, Tsukuba, 305-0801, Japan 

(Dated: February 2, 2008) 

Experimental evidence of the d-f Coulomb interaction responsible for the resonant x-ray scattering 
(RXS) from antiferroquadrupolar order in DyB2C2 is presented. The energy dependences of the 
RXS intensity with polarization analysis are analyzed by considering the interference between the 
resonances of dipolar (-E1) and quadrupolar (E2) transition processes. It is found that the structure 
factors for the El and E2 processes have the same sign for a-ir' but the opposite sign for cr- 
a' channel. This result, when compared with the calculated structure factors, means that the 
quadrupolar moments of the 4/ and 5d electrons have opposite signs to each other. Interference 
between nonresonant Thomson scattering from the lattice distortion and the resonant scatterings is 
also studied and a direct coupling between 4/ and the lattice is concluded. 

PACS numbers: 75.25.+Z, 61.10.Eq, 71.20.Eh, 75.40.Cx 



Resonant x-ray scattering (RXS) has been utilized ex- 
tensively as a powerful tool to investigate periodic align- 
ments of local anisotropics caused by, for example, chem- 
ical bondings and magnetic orderings P, . Element and 
electronic-shell selectivity provides valuable information 
on the anisotropy of a specific ion in solids. When the 
energy of an incident photon is tuned to an absorption 
edge of the element in study, a core electron is promoted 
into one of the unoccupied shells to form an intermedi- 
ate state. When it decays back to the initial state, a 
secondary photon is emitted. The cross-section of this 
resonant process carries information on the anisotropy of 
the unoccupied shell. Since it was found in manganese 
oxides that RXS can also be apowerful tool to study 
orderings of electronic orbitals [3j , the method has also 
been applied to antiferroquadrupolar (AFQ) orderings in 
/ electron systems; e.g., in DyB2C2 and in CeBg 

Observations of the orbital orderings in Mn oxides by 
RXS have been performed to date mainly using the E\ 
resonance of the K edge, which measures the anisotropy 
of the 4p state. Different theories have been proposed 
to explain the mechanism of the observation. One the- 
ory states that the 4p state reflects the anisotropy of 
the 3d state through Coulomb interaction |^. How- 
ever, since the orbital orderings in transition-metal oxides 
generally accompany Jahn- Teller lattice distortions, the 
4p states are affected also by the displacements of oxy- 
gen ions at the octahedral sites. Another theory states 
that this effect is much stronger than the Coulomb in- 
teraction because of the extended wavefunction of the 
Ap state 0- Since both theories give the same polariza- 
tion and azimuthal-angle dependences, it has been diffi- 
cult to distinguish the two interpretations. Recently, it 
was shown that the latter mechanism is dominant in Mn 
oxides from the experimental results on thin films 
In 4/ electron systems, on the other hand, the relation 



among 4/ orbital (quadrupolar moment), 5d orbital, and 
lattice has not yet been clarified. Three types of inter- 
actions should be considered: coupling between 4/ and 
lattice, bd and lattice, and d-f Coulomb interaction. In 
the present paper we focus on the RXS from AFQ order 
in DyB2C2. A theoretical estimation of the RXS inten- 
sity claims that the effect of the lattice distortion of B 
and C atoms on the Sd state is much larger than the d-f 
Coulomb interaction 

In DyB2C2 the AFQ order takes place at Tq = 24.7 
K, followed by an antiferromagnetic (AFM) order at 
Tjv — 15.3 K The magnetic structure below 

determined in Ref. is illustrated in Fig. The mo- 
ments lie in the basal c plane and those at the corner 
site and the face center site are coupled antiferromag- 
netically but canted from the (110) directions. The 
adjacent moments along the c axis are almost perpen- 
dicular with each other. This structure is characterized 
by four q-vectors: qi = {l 0), q2 = {l 0.5), q3=(0 0), 
and g4=(0 0.5). This unusual magnetic structure can 
be understood only by considering the underlying AFQ 
order, which is illustrated by the ellipses in Fig. ^ 

RXS experiment on DyB2C2 to observe the AFQ order 
was first performed by Hirota et al. and by Tanaka et al. 
using the Lm edge of Dy 0. Later, we made a detailed 
analysis on the observed resonant reflections using a the- 
ory developed by Blume [llj]. The calculated scattering 
amplitudes for each q-vector, polarization channel {a-a' 
or (t-tt'), and transition process {El or E2) assuming the 
AFQ and AFM structures shown in Fig. are mostly 
consistent with the experimental results. In short, the 
AFQ order is observed at 94 and 92, and a staggered 
lattice distortion at 92- 

However, in the previous works of RXS in DyB2C2, 
the (0 half-integer) reflections for the cr-a' channel in 
the E2 process are hidden in the tail of the El reso- 
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FIG. 1: Top view of the crystal structure of DyB2C2. Arrows 
indicate the magnetic moments in the AFM phase. Ellipses 
indicate the schematic views of pancake-like charge distribu- 
tions of 4/ electrons below Tq. Their flat planes are perpen- 
dicular to the magnetic moments. The B-C layers are located 
at z = 0.5. 
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nance, while those for the ct-tt' channel are well resolved. 
Since the calculated structure factors for the two chan- 
nels have the same order of magnitudes, the disappear- 
ance of the E2 resonance in a-a' has been a mystery. 
Recently, Tanaka et al. made a systematic RXS study of 
the (0 half-integer) reflections and analyzed the data 
using the structure factor of the atomic tensors. They 
ascribed the disappearance to the idea that the rank 2 
(quadrupole) and rank 4 (hexadecapole) tensors cancel 
out in the structure factor of the a-a' channel [l^ . 

In the present paper, we report our new results of the 
RXS for the (0 2.5)(e 94) and (3 1.5) (e 92) reflec- 
tions and show that the above problem can be solved by 
considering the interference between the El and E2 res- 
onances; the 4/-5(i Coulomb interaction will be revealed. 
Further, in the (3 1.5) reflection, the interference of the 
resonances with the Thomson scattering from the lat- 
tice distortion allows determination of the relationship 
between the orbitals and the lattice distortion; it will be 
shown that the lattice is coupled with the 4/ orbital, and 
not with the 5d orbital. 

RXS experiments were performed at the beam line 
16A2 of the Photon Factory, KEK. The incident energy 
was tuned to the Lni absorption edge of Dy and the 
polarization of the diffracted beam was analyzed using 
the PG (0 6) reflection. Flat (0 1) and (2 1) sur- 
faces of the DyB2C2 samples were prepared for the two 
experiments on the (0 2.5) and (3 1.5) reflections, 
respectively. The azimuthal angle ^' is defined to be zero 
when the b axis is along k -\- k', the vector sum of the 
incident and scattered x rays. 

Figure|21shows the energy dependences of the intensity 
of the (0 2.5) reflection in the AFQ phase for a-a' and 
ct-tt' channels. Solid lines are the fits with the structure 
factor 
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FIG. 2: Incident energy dependences of the intensity of the 
(0 2.5) reflection in the AFQ phase for the a-a' channel 
at ^' = 45° and a-n' channel at ^' = 0°. The data are cor- 
rected for the absorption. Solid lines are the fits considering 
the interference between El and E2 resonances. Dotted and 
dashed lines are the contributions from El and E2 processes, 
respectively, without the interference. 

where I-Faa'P were compared with the data. The param- 
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eters obtained are the foUowings: F. 
F^^P = 0.0396, F^^P = 0.00738, F^^P 
Ae2 = 7.782 keV, A^i = 7.790 keV, r£;2 = 2.75 eV, and 
r^i = 4.39 eV. The mixing parameter r, introduced in 
Ref. 0, was fixed at unity, which means that the El and 
E2 resonances interfere coherently. It should be noted 



that the data can be explained with finite F^^, of the 
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same order of magnitude as F„„i ' if the sign of F 

I E\\ 

opposite to that of F^^, ; this is due to the interference 
effect. 

The calculated structure factors of the even rank 
atomic tensor for the (00 half- integer) reflections, using 
the formalism in Refs.^3 and Q, are written as 



^sin2*sin20(r2 ysrf 
A cos 2^- sin 9 sin 2(j){T!f'^)^d 
for the El process and 
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for the E2 process, where sin^ 6* = 0.315 for (0 2.5). 
A and B are the positive constant factors in F^^^^ and 



respectively. Note that the El and E2 resonances 
have the same azimuthal- angle dependence for both po- 
larization channels. {Tq^"^) represents the matrix element 
of the gth component of the atomic tensor of rank K in 
the local ionic coordinates, where the principal axis of 
the quadrupolar moment is taken as the x axis and the c 
axis of the crystal as the z axis. The x axis coincides with 
the direction of the magnetic moment in the AFM phase, 
which is canted from the (10 0) directions by an angle 4>. 
Since the ion is in the local symmetry of 2 /m in the AFQ 
phase [11 m, (r,(^^) with K =even has only the com- 



ponents q = ±2 and the relation (Tj 



holds. 



Both {T^'^^) and {T^'^'^) are real and {T^'^'') represents the 
(x^ — 2;^)-type quadrupolc moment and (T2^'') represents 
the {x"^ — ~ ~ y^)?'^}-type hexadecapole moment. 
It should be noted that {Tq^^) in the structure factors 
of El and E2 processes represent atomic tensors of the 
Sd and 4/ shells, respectively, which is indicated by the 
sufBx. 



If the sign of {T^'')^^ is opposite to that of (J"2^'')4/, 
we have opposite signs for F^^^^ and F^^'^P and the same 
signs for F^^P and F^^P . Fj^^P and Fj^^P have the 
same order of magnitudes if we ignore the rank 4 ten- 
sors. This is consistent with the fitting result of the data 
and can explain the obscure and distinct features of the 
E2 resonance in the a-a' and ct-tt' channels, respectively. 
The present analysis can also explain the experimental 
results in Ref. 12. Indeed, the data for a-a' can actuall; 
be fitted with vanishing E2 term as claimed in Ref 
if we use larger Tei for a-a' . However, the rank 4 term 
in F^^^^ is not necessary to explain the result. This is 
justified by the azimuthal-angle dependence of the E2 
resonance for the (3 1.5) reflection, where only the fea- 
ture of the rank 2 tensor is observed with no indication 
of the rank 4 tensor; this point will be presented in detail 
in a forthcoming paper. 

The above analysis also applies to the RXS of the (3 
1.5) reflection, where a large nonresonant Thomson scat- 
tering from the lattice distortion also interferes. Concern- 
ing the lattice distortion below Tq, Adachi et al. recently 
presented firm evidence that B and C atoms are shifted 
along the c axis Here, two cases can be considered. 
In the first case the BC-1 parallelogram in Fig.^is shifted 
up (down) and the BC-2 is shifted down (up) at z =0.5 
(1.5), and in the second case vice versa. The structure 
factor i^iat for the two cases are the same in magnitude 
but opposite in sign. Although the sign cannot be de- 
termined by normal Thomson scattering, analysis of the 
interference with the resonance makes it possible as de- 
scribed next. 

Figure 01 shows the incident energy dependence of the 
intensity of the (3 1.5) reflection. The nonresonant 
region far below the absorption edge do not exhibit 
azimuthal-angle dependence, while the intensities at en- 
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FIG. 3: Incident energy dependences of the intensity of the (3 
1.5) reflection in the AFQ phase at * = 180° and * = 270°. 
The data are corrected for the absorption, although it is not 
perfect. The intensities well below and well above the edge 
should be equal if the correction is perfect. 



ergies of -El and E2 resonances vary with the azimuthal 
angle, indicating that there are indeed resonances. 

We have calculated the structure factors of the rank 2 
atomic tensor for the (3 1.5) reflection and have found 
that the sign of {T2 )5d must be opposite to that of 

{T!^^)if to explain the result in Fig. 13 which is con- 
sistent with the (0 2.5) reflection. Further, in order to 
explain the resonances of E2 and El processes in the a- 
a' channel, where the Thomson scattering also interferes, 

f2l 

the sign of Fjat must be opposite to that of (T2 )4/ be- 
cause the coefficient of {T!^^)if in F^^fP is negative. This 
means that the first case of the shift of the B-C parallelo- 
grams is realized, and not the second case. A simulation 
of the interference among nonresonant Thomson scatter- 
ing, E2 resonance, and El resonance is demonstrated in 
Fig. 01 Fcaic represents the coefficient of {T^'^) in the 
calculated structure factors. 

f 21 

It must be noticed that the matrix element of (T2 ) 
involves the reduced matrix element, which contains in- 
formation of the intermediate state and can be positive or 
negative depending on the ion species. When the charge 
distribution of 4/ electrons of a Dy ion is extended along 
the y direction of the local coordinates, the quadrupo- 
lar moment O22 {— Jx — Jy) is positive. The matrix 
(2'] 

element (T2 )if then becomes negative because the re- 
duced matrix element of Dy'^^ for rank 2 is negative ^3 ■ 

f 2) 

Non-zero (Tj )zd value is caused by the lifting of the 
degeneracy of yz and zx type 5d orbitals in the tetrago- 
nal environment. The local lattice distortion of B and C 
atoms in the present case will favor the yz type orbital, 
which extends along the y direction and has less mixing 
with the 2p orbitals of B and C than the zx type orbital. 



4 



800 
I 600 
-E 400 

■^L^ 200 

eg 

uf 50 




1 1 1 1 1 


I I 
[\ 1'=270° 


— u/_-|Qn° 

\ 




1 1 1 1 1 


1 1 1 1 1 1 1 


1 1 1 1 1 

— G-Jt' 


1 1 1 1 1 1 1 


1 1 1 1 1 


1'=180° : 

.yy^i , , , , 1 = 


7.75 


7.80 7.85 



Energy [keV] 



FIG. 4: Simulation of the energy dependence of the intensity 
of the (3 1.5) reflection at * = 180° and * = 270°. The 
parameters are taken to be flat = 20, F'^^^ = Fjf^^' x (-0.1), 

F'-"^^ = -P^al'' X 0.033, Vei = Tea = 6 eV, Ae2 = 7.78 keV, 
and A_Ei — 7.79 keV. The rank 2 structure factors for the 
a-n' channel vanishes at 'I' = 270° . 



However, the calculation of (T2 )5d for the yz orbital, 

following the formalism described in Ref. fl^ gives a neg- 

ative value, whereas (Tj )5d becomes positive for the zx 

orbital. Then, the result of the analysis that the signs of 

(2) (2) 
{T2 )if and (Tj )zd are opposite indicates that the zx 

orbital is favored in the state. 

The present analysis assumes localized atomic states 
both for 4/ and Sd states with fixed energy levels in the 
framework of the idealized scattering length |l3l |. We 
have also assumed that Fiat, and F'^^^ interfere 

coherently. Actually, a resonance is composed of many 
intermediate states with slightly different energies. Al- 
though the consistent explanation of the data supports 
these assumptions, it should be examined by more real- 
istic theoretical calculations that take into account com- 
plex energy levels of the intermediate state . 

In conclusion, the charge distributions of the 4/ and bd 
electrons of Dy are extended along the directions that are 
orthogonal with each other because of the d-f Coulomb 
interaction. This result is consistent with the recent re- 
examination of the RXS intensity . The relation with 
the lattice distortion is such that the B-C atoms in the di- 
rection where the 4/ charge distribution is extended move 
away from Dy, while those in the direction where the bd 
charge distribution is extended move close to Dy. Al- 
though the 5d orbital seems to be strongly coupled with 
the lattice, the present result indicates that the 4/ or- 
bital is much more strongly coupled with the lattice. The 
resultant lattice distortion is such that the quadrupole- 
strain interaction of the 4/ electrons gains energy, over- 
whelming the coupling between bd and the lattice. This 
may be because the yz and zx orbitals in the 5c? state 
are almost empty and they have little effect on the total 



energy when the lattice is distorted. Energy band calcu- 
lation of LaB2C2 shows that the conduction band of the 
La 5d state consists of the xy orbital [l7]. It is expected 
that a small amount of zx component is induced in the 
AFQ phase by the d-f Coulomb interaction, giving rise 
to the El resonance. 
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